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Quantum mechanics imposes 'monogamy' constraints on the sharing of entanglement. We show that, despite 
these limitations, entanglement can be fully 'promiscuous', i.e. simultaneously present in unlimited two-body 
and many-body forms in states living in an infinite-dimensional Hilbert space. Monogamy just bounds the 
divergence rate of the various entanglement contributions. This is demonstrated in simple families of A'-mode 
(A' > 4) Gaussian states of light fields or atomic ensembles, which therefore enable infinitely more freedom 
in the distribution of information, as opposed to systems of individual qubits. Such a finding is of importance 
for the quantification, understanding and potential exploitation of shared quantum correlations in continuous 
variable systems. We discuss how promiscuity gradually arises when considering simple families of discrete 
variable states, with increasing Hilbert space dimension towards the continuous variable limit. Such models 
are somehow analogous to Gaussian states with asymptotically diverging, but finite squeezing. In this respect, 
we find that non-Gaussian states (which in general are more entangled than Gaussian states), exhibit also the 
interesting feature that their entanglement is more shareable: in the non-Gaussian multipartite arena, unlimited 
promiscuity can be already achieved among three entangled parties, while this is impossible for Gaussian, even 
infinitely squeezed states. 

PACS numbers: 03.67.Mn, 03.65.Ud 



I. INTRODUCTION 

Entanglement is a core concept in quantum mechanics fl'], 
and a key resource for quantum communication and infor- 
mation processing 0]. Unlike classical coiTelations, entan- 
glement cannot be freely distributed |3]. In a multipartite 
compound system of two-level quantum objects (qubits), if 
two subsystems are maximally entangled, they cannot share 
any residual form of quantum correlations with the other re- 
maining parties 10, Ht]. Analogous monogamy relations have 
been recently established for entanglement between canoni- 
cal conjugate variables of continuous variable (CV) systems 
10,13], like harmonic oscillators, light modes and atomic en- 
sembles, endowed with an infinite-dimensional Hilbert space. 
In the general case of a state distributed among parties 
(each owning a single qubit, or a single mode, respectively), 
the monogamy constraint on bipartite entanglement takes the 
form of the Coffman-Kundu-Wootters inequality [4J 
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where the global system is multipartitioned in subsystems S k 
{k - 1, . . .,A^), each owned by a respective party, and £ is a 
proper measure of bipartite entanglement. The left-hand side 
of inequality ([T) quantifies the bipartite entanglement between 
a probe subsystem S , and the remaining subsystems taken as a 
whole. The right-hand side quantifies the total bipartite entan- 
glement between S and each one of the other subsystems S j^,! 
in the respective reduced states. The non-negative difference 
between these two entanglements, minimized over all choices 
of the probe subsystem, is referred to as the residual multi- 
partite entanglement. It quantifies the purely quantum corre- 



lations that are not encoded in pairwise form, so it includes all 
manifestations of genuine /T-partite entanglement, involving 
K subsystems at a time, with 2 < K < N. In the simplest non- 
trivial instance of = 3, the residual entanglement has the 
(physical and mathematical) meaning of the genuine tripartite 
entanglement shared by the three subsystems 14 6, 8]. In the 
general case A^ > 3, a natural way to discriminate between 
all those entanglement contributions has been very recently 
advanced in Ref. |9], and proven successful when address- 
ing entanglement shared by harmonic systems under complete 
permutation invariance. The study of entanglement sharing 
and monogamy constraints thus offers a natural framework to 
interpret and quantify entanglement in multipartite quantum 
systems 

From an operational perspective, qubits are the main logi- 
cal units for standard realizations of quantum information (QI) 
protocols [2]. Also CV Gaussian entangled resources have 
been proven useful for all known implementations of QI pro- 
cessing fTo*!, including quantum computation 1 11 ], sometimes 
outperforming more traditional qubit-based approaches as in 
the case of teleportation lfl2ll . It is therefore important to un- 
derstand if special features of entanglement appear in states 
of infinite Hilbert spaces, which are unparalleled in the corre- 
sponding states of qubits. Such findings may lead to new ways 
of manipulating QI in the CV setting ||13|1 . and may in general 
contribute to a deeper and more complete understanding of 
entanglement in complex systems. 

In this paper, we address this motivation on a clear-cut 
physical ground, aiming in particular to show whether the 
unboundedness of the mean energy characterizing CV states 
enables a qualitatively richer structure for distributed quan- 
tum correlations. We prove that multimode Gaussian states 
exist, that can possess simultaneously arbitrarily large pair- 



wise bipartite entanglement between some pairs of modes an 
arbitrarily large genuine multipartite entanglement among a. 
modes. In particular, we focus on a four-mode family of Gaus 
sian states which are producible with standard optical means 
the achievable amount of entanglement being technologicall 
limited only by the attainable degree of squeezing. Thes 
states asymptotically reach the form of two perfectly entail 
gled Einstein-Podolsky-Rosen (EPR) pairs that can moreove 
be arbitrarily intercorrelated quantumly, as well as share a 
asymptotically diverging genuine four-partite entanglemer 
among all modes. We show that this promiscuity is fully com 
patible with the monogamy inequality ([TJ for CV entangle 
ment as established in Refs. |6, 7], and even with the stronge 
monogamy constraints on distributed bipartite and multipai 
tite quantum correlations recently put forward in Ref. |9 
This original feature of entanglement sheds new light on th 
actual extent to which the basic laws of quantum mechar 
ics curtail the distribution of information; on a more applies 
tive ground, it may serve as a prelude to implementations c 
QI processing in a infinite-dimensional scenario that cannc 
be achieved with single qubit resources. While promiscu 
ity might seem an unexpected property of quantum correla 
tions, we will show that it indeed arises naturally and gradu 
ally with increasing dimension of the Hilbert space associate 
to the entanglement-sharing parties. In the limit of continuou 
variable systems, non-Gaussian states will be shown to be i 
general sensibly more powerful than Gaussian states from th 
point of view of (promiscuous) entanglement sharing. 



II. ENTANGLEMENT SHARING AMONG GAUSSIAN 
MODES 

Entanglement in C V systems is encoded in the form of EPR 
correlations [TT]. Let us consider the motional degrees of 
freedom of two particles, or the quadratures of a two-mode 
radiation field, where mode k - i, j is described by the lad- 
der operators a^, a'j^ satisfying the bosonic commutation rela- 
tion [at, a]l - 1. An arbitrarily increasing degree of entan- 
glement 1I15I1 (and of mean energy) can be encoded in a two- 
mode squeezed state = Uij(r) (|0),® |0)j) with increas- 
ing squeezing factor r e R, where the (phase-free) two-mode 
squeezing operator is lfT6ll Uij(r) = exp |^^(fl^!aj - 0,0;) j, and 
|0)t denotes the vacuum state in the Fock space of mode k. In 
the limit of infinite squeezing (r — > 00), the state approaches 
the ideal (unnormalizable, infinite-energy and unphysical) 
EPR state (T^, simultaneous eigenstate of total momentum 
and relative position of the two subsystems, which thus share 
infinite entanglement. The entanglement of squeezed states 
and, more generally, Gaussian states, has been intensively 
studied in recent times [17]. Gaussian states of modes 
are completely described in phase space by the 2N x 2N 
real symmetric covariance matrix (CM) cr of the second mo- 
ments J (XiXj + XjXi^ of the canonical bosonic operators 

Xk = qk ^ h + a[,XN+k = Pk ^ ak - a[(k = l...A^). 
In this representation the unitary two-mode squeezing opera- 
tor amounts to a symplectic matrix S ij{r) = (^J ^' j © ^ 
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FIG. 1: (Color online) (a) Preparation of the four-mode Gaussian 
states y of Eq. lO via two-mode squeezers (boxes) applied on vac- 
uum modes (yellow beams), (b) Structure of bipartite entanglement. 



(where = cosh r, Sr - sinh r, which acts by congruence on 
the CM, o-h^ So-S^. 

Entanglement sharing has been addressed for Gaussian 
states in Refs. by adapting and extending the original 

CofFman-Kundu-Wootters analysis [4] (see also |18] where a 
different analysis is performed, which nonetheless has impli- 
cations for the sharing of correlations in multimode Gaussian 
states). In the most basic multipartite CV setting, namely 
that of three-mode Gaussian states, a partial "promiscu- 
ity" of entanglement can be achieved. Permutation-invariant 
states exist which are the Gaussian simultaneous analogues of 
Greenberger-Horne-Zeilinger (GHZ) and W states of qubits 
||8l[l9t]. They possess unlimited tripartite entanglement (with 
increasing squeezing) and nonzero, accordingly increasing bi- 
partite entanglement which nevertheless stays finite even for 
infinite squeezing |6]. We will now show that in CV systems 
with more than three modes, entanglement can be distributed 
in an infinitely promiscuous way. To illustrate the existence 
of such phenomenon, we consider the simplest nontrivial in- 
stance of a family of four-mode Gaussian states, endowed 
with a partial symmetry under mode exchange. 



III. FOUR-MODE STATES: STRUCTURAL AND 
ENTANGLEMENT PROPERTIES 

We start with an uncorrected state of four modes, each 
one initially in the vacuum of the respective Fock space, 
whose corresponding CM is the identity. We apply a two- 
mode squeezing transformation with squeezing s to modes 
2 and 3, then two further two-mode squeezing transforma- 
tions (redistributing the initial pairwise entanglement among 
all modes) with squeezing a to the pairs of modes 1,2 and 3,4 
[see Fig.[Tla)]. For any value of the parameters s and a, the 
output is a pure four-mode Gaussian state with CM 7, 
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A state of this form is invariant under the double exchange of 
modes 1 <-> 4 and 2 «-> 3, as S ij = S jj and operations on 
disjoint pairs of modes commute. 

In a pure four-mode Gaussian state and in its reductions, bi- 
partite entanglement is equivalent to negativity of the partially 
transposed CM, obtained by reversing time in the subspace 
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of any chosen single subsystem 11201 12111 . This inseparability 
criterion is readily verified for the family of states in Eq. ^ 
yielding that, for all nonzero values of the squeezings s and 
a, y is entangled with respect to any global bipartition of the 
modes. The state is thus said to he fully inseparable 1^, i.e. 
it contains genuine four-partite entanglement. Following pre- 
vious studies on CV entanglement sharing we use the "con- 
tangle" T |6] to quantify bipartite entanglement, an entangle- 
ment monotone under Gaussian local operations and classical 
communication. It is defined for pure states as the squared 
logarithmic negativity ll22ll (which quantifies how much the 
partially transposed state fails to be positive) and extended to 
mixed states via Gaussian convex roof ll23i l24ll . i.e. as the 
minimum of the average pure-state entanglement over all de- 
compositions of the mixed state in ensembles of pure Gaus- 
sian states. If (Tiy is the CM of a (generally mixed) bipartite 
Gaussian state where subsystem / comprises one mode only, 
then the contangle r can be computed as 

Tio-iij) = t{<t°^I) = g[mlj\, g[x\ = arcsinh^[ ^x - 1]. (3) 

Here cr°'j" corresponds to a pure Gaussian state, and m,|j = 

m(cr"y') - -^Det cr"'" - ^Det cr^', with (r"IJ^ being the re- 
duced CM of subsystem obtained by tracing over the de- 
grees of freedom of subsystem j (/). The CM cr°^' denotes the 

pure bipartite Gaussian state which minimizes m((r'^y) among 
all pure-state CMs cr^^. such that crjj^. < (T,|j. If (riy is a 
pure state, then cr'!y' - criy, while for a mixed Gaussian state 
Eq. OJ is mathematically equivalent to constructing the Gaus- 
sian convex roof. For a separable state m(tr°^') = 1. The 
contangle t is completely equivalent to the Gaussian entan- 
glement of formation ll23ll . which quantifies the cost of creat- 
ing a given mixed Gaussian state out of an ensemble of pure, 
entangled Gaussian states. 



A. Structure of bipartite entanglement 



S ij(a)S Jj(a) of modes / and j generated with the same squeez- 
ing a. In fact, the two-mode mixed state yip (and, equiva- 
lently, 7314) serves as a proper resource for CV teleportation 
lfT2ll . realizing a perfect transfer (unit fidelity |25]) in the limit 
of infinite squeezing a. The four-mode state 7 reproduces 
thus (in the regime of very high a) the entanglement con- 
tent of two EPR-like pairs (1,2 and 3,4). Remarkably, there 
is an additional, independent entanglement between the two 
pairs, given by t(7(12)|(34)) = [the original entanglement in 
the two-mode squeezed state S2,^{s)S^t(^s) after the first con- 
struction step, see Fig.lTfa)], which can be itself increased ar- 
bitrarily with increasing s. This peculiar distribution of bipar- 
tite entanglement [see Fig.[TJb)] is a first interesting signature 
of an unmatched freedom of entanglement sharing in multi- 
mode Gaussian states as opposed for instance to states of the 
same number of qubits, where a similar situation is impossi- 
ble. Specifically, if in a pure state of four qubits the first two 
approach unit entanglement and the same holds for the last 
two, the only compatible global state reduces necessarily to a 
product state of the two singlets: no interpair entanglement is 
allowed by the monogamy constraint HB- 



B. Residual entanglement 

We can now move to a closer analysis of entanglement 
distribution and genuine multipartite quantum correlations. 
A first step is to verify whether the monogamy inequal- 
ity ([T]i is satisfied on the four-mode state 7 distributed among 
the four parties (each one owning a single mode) ll26ll . In 
fact, the problem reduces to proving that min{g[OTju.234j] - 

^['w^p]' ^['"21(134)] - ^['"^p] - ^['"213]' nonnegative. The 
first quantity always achieves the minimum yielding 



t''''(7) = t(7,|(234))-t(7ip) 



(4) 
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arcsinh < ^^ [cosh a + cosh(2s) sinh fl]^ - 1 > - 4fl 



In the four-mode state with CM y,'^& can compute the 
bipartite contangle in closed form 02411 for all pairwise re- 
duced (mixed) states of two modes i and j, described by a 
CM 7,|y. By applying the partial transpose criterion [20], 
we find that the two-mode states indexed by the partitions 
1|3, 2|4, and 1|4, are separable. For the remaining two- 
mode states we find mip - - cosh2fl, while m2|3 

is equal to -l+2':°^h^P'')':°^h^,+3cosh(2.0-4sinh^«sinh(2.,) jf ^ < 
4[cosh" a+e^'' sinh a] 

arcsinh[ Vtanh s], and to 1 (implying separability) otherwise. 
Accordingly, one can compute the pure-state entanglements 
between one probe mode and the remaining three modes. 
One finds mi|(234) = ?M4|(i23) = cosh^ a + cosh(2s) sinh^ a and 
W2|(i34) — 'M3|(i24) = sinh^ a + cosh(2s) cosh^ a. 

Concerning the structure of bipartite entanglement, the con- 
tangle in the mixed two-mode states 7ip and 7314 is Aa^, un- 
respective of the value of s. This quantity is exactly equal to 
the degree of entanglement in a pure two-mode squeezed state 



Since cosh(2.s) > 1 for 5 > 0, it follows that t''''^ > 0. The en- 
tanglement in the global Gaussian state is therefore correctly 
distributed according to the monogamy law, in such a way 
that the residual contangle t'''"* quantifies the multipartite en- 
tanglement not stored in couplewise form. Those quantum 
correlations can be either tripartite involving three of the four 
modes, and/or genuinely four-partite among all of them. Con- 
cerning the tripartite entanglement, we first observe that in the 
tripartitions 1|2|4 and 1|3|4 the tripartite entanglement is zero, 
as mode 4 is not entangled with the block of modes 1,2, and 
mode 1 is not entangled with the block of modes 3,4 

The only tripartite entanglement present, if any, is equal 
in content (due to the symmetry of the state 7) for the triparti- 
tions 1|2|3 and 2|3|4, and can be quantified by the residual con- 
tangle (a Gaussian entanglement monotone |6]) determined 
by the corresponding three-mode monogamy inequality ([1]). 
One can show (see caption of Fig. |2]i that such genuine tri- 
partite entanglement T'"(7ip|3) is bounded from above by the 
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FIG. 2: (Color online) Upper bound t''°""''(7ip|3), Eq. lO, on the tri- 
partite entanglement between modes 1, 2 and 3 (and equivalently 2, 
3, and 4) of the four-mode Gaussian state defined by Eq. lO, plotted 
as a function of the squeezing parameters s and a. Focusing on the 
tripartition 1|2|3, the bipartite contangle t(7,|(^q) (with i, j, k a permu- 
tation of 1,2,3) is bounded from above by the corresponding bipartite 
contangle T(t7^j^.j.j) in any pure, three-mode Gaussian state with CM 
°'l(]K) ^ r,io«- The state a" = S i2{a)S 2,i{t)Sl^(t)Sl^{a) of the three 
modes 1, 2 and 3, with / = ^arccosh[ ' ], satisfies this con- 
dition and from Eq. l|3j we have Tiy^ji^^) < g[(raj[°™)'')^]. This leads 
to Eq. (|5](, where the quantity g[(in''r^""^^)^} — g[m\^^ — g[m^^j] with 

ra*°""^ = sinh^ a + mjjj™'' cosh" a is not included in the minimization, 
being always larger than the other terms. 



quantity 

t'"'""''(7i|2|3) = min{g[(m?-«^)2]-^[m22], g[(m'^^,p^]-g[ml,]}, 

(5) 

with = l+^e'^h^^tanh^^ hound ^ ^Qgj^2 ^ ^ ^^bouml ^^^^^2 ^ 

3|(12) 1-sech^fltanh^ .5 11(23) 3|(12) 

The upper bound t^"""'' {y [\2n) is always nonnegative (as a 
consequence of monogamy |6]), is decreasing with increas- 
ing squeezing a, and vanishes in the limit a — > oo, as shown 
in Fig. (|2]l. Therefore, in the regime of increasingly high a, 
eventually approaching infinity, any form of tripartite entan- 
glement among any three modes in the state 7 is negligible 
(exactly vanishing in the limit). As a crucial consequence, in 
that regime the residual entanglement t"'(7) determined by 
Eq. (|4]l is all stored in four-mode quantum coiTelations and 
quantifies the genuine four-partite entanglement. 

We finally observe that VHy) is an increasing function of 
a for any value of s (see Fig. |3]l, and it diverges in the limit 
fl — > 00. This proves that the class of pure four-mode Gaus- 
sian states with CM 7 given by Eq. ^ exhibits genuine four- 
partite entanglement which grows unboundedly with increas- 
ing squeezing a and, simultaneously, possesses pairwise bi- 
partite entanglement in the mixed two-mode reduced states 
of modes 1,2 and 3,4, that increases unboundedly as well 
with with increasing a f\3\. Moreover, as previously shown, 
the two pairs can themselves be arbitrarily entangled with 
each other with increasing squeezing s. Notice that usual 
monogamy inequalities do not typically constrain "hybrid" 
entanglement distributions such as 2 x 2 versus 1 x 1 or 1 x 3, 
even though it is reasonable to expect that some limitations 
exist also when entanglement is shared over such partitions. 



FIG. 3: (Color online) Residual multipartite entanglement T'"{y) 
[Eq. Q], which in the regime of large squeezing a is completely 
distributed in the form of genuine four-partite quantum correlations. 



C. Strong monogamy and genuine four-partite entanglement 

It would be desirable to have a clear distinction between tri- 
partite and four-partite quantum correlations also in the range 
in which the former are not negligible, i.e., for small a. A 
simple idea comes from Ref. |9], which adapted to our set- 
ting states the following. The minimum residual entangle- 
ment Eq. (HJi between one probe mode and the remaining 
three, is equal to the total three-mode entanglement involv- 
ing the probe mode and any two of the others, plus the gen- 
uine four-partite entanglement shared among all modes. Such 
a construction, which provides an additional decomposition 
for the difference between LHS and RHS of Ineq. ([T]), leads 
in general to postulate that a stronger monogamy constraint 
acts on entanglement shared by parties, which imposes a 
trade-off on both the bipartite and all forms of genuine multi- 
partite entanglement on the same ground. Existence of such 
a constraint, which directly generalizes the Coffman-Kundu- 
Wootters (A^ = 3) case |4], has been indeed proven for ar- 
bitrary Gaussian states endowed with permutation invariance 
(i.e., for the generalized A^-mode Gaussian GHZ/W states and 
their mixed reductions), and a bona fide quantification of the 
genuine A^-party entanglement (computed again by means of 
the contangle) has been obtained on such states 

We will now show that strong monogamy also holds in 
the partially symmetric four-mode Gaussian states of Eq. 
which is an important indication of the generality of the ap- 
proach presented in [9]. Indeed, strong monogamy imme- 
diately follows from Eqs. (I4l5l l, as T''''Hy) - t"''(7i|2|3) > 
T'-'^'iy) > t''""'"'(7i|2|3) = T''*''""''(ri|2|3|4) > 0. The quantity 
^Abound ^^^^^^^^ jg ^ /ower bound to the genuine four-partite en- 
tanglement shared by the four modes in the state 7, for any a 
and s. Its functional dependence on these squeezing parame- 
ters is very similar to that of t'''^\'y) (one needs only to sub- 
tract the surface in Fig. |2l which extends over a very narrow 
scale, from the surface in Fig. [3]). r"°*(7) represents actually 
an upper bound for the four-partite entanglement, which gets 
asymptotically tight for a » 1: in this limit, upper and lower 
bounds coincide and as already remarked, the residual entan- 
glement is all in four-partite form. This brief analysis is very 
relevant on one hand for the purposes of Ref. |9] as it embod- 
ies the first proof of strong monogamy beyond strong sym- 
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metry requirements; while on the other hand, in the scope of 
the present paper, it enables us to regard the considerations of 
the previous subsection as mathematically valid over all the 
parameter space, and not only in the limit of high squeezing 
a. 



IV. DISCUSSION: WHERE DOES PROMISCUITY COME 
FROM? 

A. Some operational considerations 

From a practical point of view, two-mode squeezing trans- 
formations are basic tools in the domain of quantum optics 
lfl6ll (they occur e.g. in parametric down-conversions), and 
the amount of producible squeezing in experiments is con- 
stantly improving [27]. Only technological, no a priori limi- 
tations need to be overcome to increase a and/or s to the point 
of engineering excellent approximations to the demonstrated 
promiscuous entanglement structure in multimode states of 
light and atoms (see also [i2&]). To make an explicit exam- 
ple, already with realistic squeezing degrees like s = 1 and 
a = 1.5 (corresponding to ~ 3 dB and 10 dB, respectively), 
one has a bipartite entanglement of t(7i|2) = ^(7314) - 9 ebits 
(corresponding to a Gaussian entanglement of formation [, 23,] 
of ~ 3.3 ebits), coexisting with a residual multipartite entan- 
glement of t"'*(7) ^5.5 ebits, of which the tripartite portion 
is at most t''"""'' (y H2\i) - 0-45 ebits. This means that one can 
simultaneously extract at least 3 qubit singlets from each pair 
of modes (1,2) and (3, 4), and more than a single copy of gen- 
uinely four-qubit entangled states (like cluster states). Albeit 
with imperfect efficiency, this entanglement transfer can be re- 
alized by means of Jaynes-Cummings interactions [29], rep- 
resenting a key step for a reliable physical interface between 
fields and qubits in a distributed QI processing network. 



B. Emergence of promiscuity from discrete to continuous 
variable systems, and beyond Gaussian states 

By constructing a simple and feasible example we have 
shown that, when the quantum correlations arise among de- 
grees of freedom spanning an infinite-dimensional space of 
states (characterized by unbounded mean energy), an accord- 
ingly infinite freedom is tolerated for QI to be doled out. The 
construction presented here can be straightforwardly extended 
to investigate the increasingly richer structure of entanglement 
sharing in A^-mode (N even) Gaussian states via additional 
pairs of two-mode squeezing operations which further redis- 
tribute entanglement. Notice once more that the promiscuous 
entanglement distribution happens with no violation of the a 
priori monogamy constraints that retain their general validity 
in quantum mechanics. The motivation is the following. In 
the CV scenario, naively speaking, if one writes from Ineq. ([T]l 
that the entanglement between one mode and the rest is equal 
to the sum of bipartite entanglements plus the residual multi- 
partite entanglement, one can in principle approach the limit- 
ing 'fully promiscuous' case "00 = 00-1- 00". The only restric- 



tion imposed by monogamy is to bound the divergence rates 
of the individual entanglement contributions as the squeezing 
parameters are increased. Within the restricted Hilbert space 
of four or more qubits, instead, an analogous entanglement 
structure between the single qubits is strictly forbidden (one 
can either have "1 = 1 -1- 0" or "1 = H- !")• 

Our results open original perspectives for the understanding 
and characterization of entanglement in multiparticle systems. 
Gaussian states with finite squeezing (finite mean energy) are 
somehow analogous to discrete systems with an effective di- 
mension related to the squeezing degree. As the promiscuous 
entanglement sharing arises in Gaussian states by asymptoti- 
cally increasing the squeezing to infinity, it is natural to expect 
that dimension-dependent families of states will exhibit an 
entanglement structure gradually more promiscuous with in- 
creasing Hilbert space dimension towards the C V limit. Let us 
further clarify this point. In the bipartite, pure-state instance, 
a two-mode squeezed Gaussian state |iAj^('")) with squeezing 
r = (logt/)/2 has exactly the same entanglement E = logii 
(quantified by the logarithmic negativity [22]) as a maximally 
entangled state of two qut/its, lO^) - d^^'^ Y^ti I") ^ 
Thus, for any finite squeezing, Gaussian entanglement is not, 
in a sense, taking advantage of the full infinite-dimensional 
Hilbert space, but only of a restricted, finite section of it. Only 
for r — > 00, and accordingly c/ — > 00, the states \i//sqir — > 00)) 
and l^d^co) tend towards the (unphysical) EPR state, and 
the entanglement and the mean energy asymptotically diverge 
id. 

With such an analogy in mind, let us now construct the sim- 
plest conceivable family of tripartite states, in which we will 
witness promiscuity of distributed entanglement arising with 
increasing dimensionality of the Hilbert space. We define the 
pure permutationally-invariant state \^2n) g of a system 
of 3A^ qubits (where A^ > 2 is an even integer, and we label 
the state with the subscript "2A" for future convenience) as 
the tensor product of N/2 copies of the three-qubit GHZ state 
ifTgll . and of N/2 copies of the three-qubit W state [8]. In for- 
mula, 

mN)^\l/^onz)''^ ®\<J^„r^ , (6) 

where |<Ac„z> = (|000) + 1111»/ V2 and \^^} = (|001) + 1010) + 
1 100))/ V3. By grouping our 3N qubits into three parties A, B 
and C, each formed by N qubits (each qubit, in turn, taken 
from a single copy of either the GHZ or the W state), the 
state e C'' » C"' » C' of Eq. ^ becomes the ti'ipar- 
tite state of the three qut/its ABC with d = 2N. A graphi- 
cal depiction of such a construction is provided in Fig. |4] It 
is straightforward now to evaluate the entanglement proper- 
ties of l^rf). Using the tangle as an entanglement measure 
ll30ll we find that the three quc/its ABC share a genuine tri- 
partite residual entanglement T^^^^'~ = d/4, as the GHZ state 
has unit three-tangle, and the W state has zero three-tangle 
[13, Bl- Similarly, the bipartite entanglement between any two 
quc/its is T'j*'* - t^'*" = r^'*' - d/9, as the two-qubit reduc- 
tions of W states have a tangle of 4/9, while the two-qubit 
reductions of GHZ states are separable. Monogamy is clearly 
satisfied (the tangle between one party and the other two is 
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FIG. 4: (Color online) A tripartite pure state of three qurfits {d = 2N 
with A' even) A, B, and C, is constructed as the tensor product of N/2 
copies of a three-qubit GHZ state, and N/2 copies of a three-qubit 
W state. In the limit d oo, ABC constitute a tripartite CV system 
whose (non-Gaussian) state exhibits full promiscuity, i.e. unlimited 
genuine tripartite entanglement and unlimited bipartite entanglement 
between any two parties. 



rj-A\(BC) _ ij 1^/26), and one can see that both - tripartite 
and every couplewise - forms of entanglement are linearly 
increasing with the dimension d, and are so monotonically 
increasing functions of each other Bipartite and genuine mul- 
tipartite entanglement can grow unboundedly (i.e., diverging 
for d — > oo), and are mutually enhanced: this is exactly what 
the concept of promiscuity means. Only for a three-qubit sys- 
tem there is mutual exclusion between maximum tripartite en- 
tanglement (the GHZ case) and maximum two-party reduced 
entanglement (the W case). As soon as quc/its with d > 2 
are concerned, an increasingly more freedom is let for entan- 
glement distribution compatibly with monogamy. An alterna- 
tive demonstration of the full promiscuity of entanglement in 
the state of Eq. (|6]l is provided in Appendix lAlbv employing 
the 'squashed entanglement' |31] as bipartite entanglement 
monotone. 

From this simple example we learn several things. First, 
that promiscuity gradually arises with increasing dimension d 
of the Hilbert space, mediating between the prohibitive qubit 
scenario, and the libertine CV one. Second, that the analogy 
between finite-squeezing CV states and finite-dimensional 
systems can be extended to the multipartite scenario. Third, 
and perhaps more interestingly, we see that Gaussian states are 
too structurally simple as compared to general non-Gaussian 
states of CV systems, and are hence sometimes unable to 
achieve the full power of entanglement distribution in infinite- 
dimensional systems even in the limit of a diverging mean 
energy per mode. What we mean is quite straight. Taking 
the limit li ^ oo in the above qut/it picture, the state \^'d^cc} 



exhibits asymptotically diverging entanglement both in gen- 
uine tripartite form and in couplewise form between any two 
parties of the CV systems A, B and C, i.e., an unlimited, 
full promiscuous structure. We know instead that three-mode 
Gaussian states can never achieve such a distribution of entan- 
glement (the reduced two-mode entanglement stays finite 
even for infinite squeezing), and in this paper we had to resort 
to four-mode states with broken symmetry to demonstrate un- 
limited promiscuity in the Gaussian scenario. 

The state |*F^^co) is indeed a highly non-Gaussian state: a 
recently introduced quantifier of non-Gaussianity |32], which 
measures the normalized Hilbert-Schmidt distance 6 between 
a non-Gaussian state and a reference Gaussian state with the 
same first and second moments, yields llssll 



(7) 



On a scale ranging from zero to one, 6{\^d)) ^ 0.48 for any 
d > A (recall that d - A corresponds to just one copy of the 
GHZ state and one copy of the W state), and increases with 
d converging to 1/2 for d ^ oo. As a term of comparison, 
1/2 is the asymptotic non-Gaussianity of the n* single-mode 
Fock state \n) for n » 1 |32]. We have thus given the first 
evidence that monogamy of CV entanglement holds beyond 
the Gaussian case (although the state we considered is highly 
non generic, see the discussion in ll34ll ). and that, crucially, 
non-Gaussian states can achieve a strictly more promiscuous 
structure in the sharing of quantum correlations, than Gaus- 
sian states. More promiscuity means more and stronger forms 
of entanglement encoded in different multipartitions and si- 
multaneously accessible for quantum information purposes. 

One can readily provide generalizations of the above ex- 
ample to convince oneself that for an arbitrary number of 
quiiits in the limit d ^ oo, & fully promiscuous structure of 
entanglement is available (under and beyond permutation in- 
variance) for the corresponding multipartite non-Gaussian C V 
states. Permutation-invariant Gaussian states (the so-called 
Gaussian GHZ/W states |6]), instead, are always promiscu- 
ous for any A^, but never fully promiscuous, in the sense that 
the reduced bipartite and/or A'-partite entanglements (K < N) 
in the various partitions always stay finite (and increasingly 
low with higher A^) even in the limit of infinite squeezing |9]. 
And we stress that this is not a consequence of restricting to 
a single mode per party: even considering a big MA^-mode 
CV system, partitioned into A^ subsystems each comprising 
M modes, the corresponding A^-partite Gaussian entanglement 
(under permutation invariance) is exactly scale invariant, i.e., 
it is equal to the entanglement among A^ single modes (i.e. to 
the case M = 1) |9], hence can never be fully promiscuous. 
Of course, while Gaussian states cannot reach the shareability 
of non-Gaussian states, mimicking the formers via the latters 
is always possible. If one wishes to reproduce the A^-mode 
permutation-invariant Gaussian instance in terms of quc/its 
with varying d, it is sufficient e.g. for A^ = 3 to modify the 
construction of Fig.|4]by unbalancing (as a function of d) the 
fraction of GHZ versus W copies, and keeping the latter num- 
ber finite while allowing the former to diverge with d. In fact, 
as just shown, this Gaussian-like setting represents clearly not 
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the best entanglement sharing structure that states in infinite- 
dimensional Hilbert spaces may achieve. It is worth point- 
ing out, however, that in the current experimental practice, 
should a certain degree of promiscuity be required, it would 
be surely easier to produce and manipulate Gaussian states of 
four modes, with a quite high squeezing, than tensor products 
of many three-qubit states. 



C. Concluding remarks and Outlook 

Inspired by the above discussion on the origin of entangle- 
ment promiscuity, a more extended investigation into the huge 
moat of quc/its appears as the next step to pursue, in order to 
aim at developing the complete picture of entanglement shar- 
ing in many -body systems fl]. Here, we have initiated this 
program by establishing a sharp discrepancy between the two 
extrema in the ladder of Hilbert space dimensions: namely, 
entanglement of CV systems in the limit of infinite mean en- 
ergy has been proven infinitely more shareable than that of 
individual qubits. We have also given indication that such a 
transition between these two diametrically opposite cases is 
not sharp, but smeared on the whole scale of Hilbert space 
dimensionalities 2 < li < oo, as promiscuity gradually and 
smoothly arises with increasing d. 

Once a more comprehensive understanding will be avail- 
able of the distributed entanglement structure in high- 
dimensional and CV systems (also beyond Gaussian states), 
the interesting task of devising new protocols to translate such 
potential into full-power QI processing implementations can 
be addressed as well. In this respect, it is worth remark- 
ing again that non-Gaussian states, which by virtue of the 
extremality theorem llssll are systematically more entangled 
that Gaussian ones at fixed second moments in a bipartite sce- 
nario, appear also to have a richer structure when entangle- 
ment sharing is addressed. In the tripartite case, for instance, 
we have constructed an example of a fully promiscuous non- 
Gaussian state, while promiscuity of three-mode Gaussian en- 
tanglement is always partial. This provides a novel, additional 
motivation for further characterizing non-Gaussian entangle- 
ment also in the multipartite scenario, together with its possi- 
ble exploitation, e.g., to perform highly efficient teleportation 
networks t36l1 or to solve the Byzantine agreement problem 

iH. 
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APPENDIX A: NON-GAUSSIAN (FULL) VERSUS 
GAUSSIAN (PARTIAL) PROMISCUITY OF TRIPARTITE 
ENTANGLEMENT SHARING VIA THE SQUASHED 
ENTANGLEMENT 

Here we provide an alternative comparison (and contrast) 
between the structure of entanglement sharing in the permu- 
tationally invariant pure CV state obtained by taking the limit 
c/ — > oo of the three-quc/it state in Eq. and the permu- 
tationally invariant pure three-mode Gaussian states (Gaus- 
sian GHZ/W states) introduced in Ref. Even though 
more abstract, this analysis is perhaps more conclusive as 
the same entanglement measure is employed for both states. 
The squashed entanglement Esq is an interesting information- 
theoretic entanglement monotone which is, among the other 
properties, additive on tensor product states [3 1 ] . Moreover, it 
has been proven to fulfill monogamy, in the sense of Ineq. ([1]), 
for arbitrary A^-partite systems, unregardingly of the associ- 
ated Hilbert space dimensions |38]. Unfortunately, this mea- 
sure is uncomputable for almost all entangled states, even 
though it is known to be upper bounded by the entangle- 
ment cost, and lower bounded by the distillable entanglement 
(hence equal to the Von Neumann entropy of entanglement on 
pure states). Let us bound the different entanglement contri- 
butions in the state of Eq. ^ in terms of squashed entangle- 
ment. The 1x2 entanglement between one qut/it, say A, and 
the other two, is 

£ff ^^(1*?^)) = («'/4)[£ll<^c\|^,„,)) + £ff c>(|^,))] ^ 0A8d 

(Al) 

(in this case the squashed entanglement reduces on each copy 
to the computable Von Neumann entropy of one qubit). Con- 
cerning two-party entanglement, the GHZ states have all sep- 
arable two-qubit reductions; on the other hand, the bipartite 
squashed entanglement in the reduced two-qubit states of each 
W copy amounts to some quantity, say ai, which is strictly 
positive, as every entangled two-qubit state has a nonzero dis- 
tillable entanglement (PPT criterion |39]). The total bipartite 
squashed entanglement between any two quiiits in the state 
Eq. (|6]l is thus 

^sf (l*l'rf» = <^d/4 , w > . (A2) 

From the general monogamy of the squashed entanglement, 
we know that the residual three-partite entanglement is well 
defined and positive. 

Exploiting in particular the fact that squashed entanglement 
is monogamous on each W copy, E^jj-^'^^di^^}) > 2a), we im- 
mediately obtain a lower bound on the three-qut/it residual 
entanglement, 

Eff^''(\'i'd})>d/4, (A3) 

where we have used the fact that for each GHZ copy 
^'sf'^d'/'GHz)) = 1- Looking at Eqs. ( IATM3] i, we observe 
once more that all forms of (squashed) entanglement are mu- 
tually increasing functions of each other on the state of Eq. (|6]l. 



8 



and diverging in the limit d ^ oo. This definitely proves that 
the corresponding asymptotic (non-Gaussian) C V state, is def- 
initely fully promiscuous, compatibly with the monogamy of 
entanglement sharing. 

On the other hand, let us repeat the same analysis for 
the Gaussian counterparts, represented by pure permutation- 
invariant three-mode "GHZ/W Gaussian states [6] with 
asymptotically diverging squeezing (and mean energy) in any 
single mode. The 1x2 squashed entanglement, equal to the 
Von Neumann entropy of one mode, diverges with the squeez- 
ing (for r » it goes ~ 2r). The reduced two-mode squashed 



entanglement, instead, is smaller than the corresponding en- 
tanglement of formation (which equates the entanglement cost 
as its additivity is proven in this special case The latter 

increases with the squeezing, converging to the finite value 
» 0.278 in the limit r — > oo. This unambiguosly prove 
that permutation-invariant tripartite Gaussian states exhibit 
only a partial promiscuity, as the genuine tripartite entangle- 
ment (obtained by difference) diverges for infinite squeezing, 
while the reduced bipartite entanglement stays constant. It is 
also known that no three-mode Gaussian states can be more 
promiscuous than these GHZ/W states ||^. 
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